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Abstract. - We study thermoelectric properties of inhomogeneous nanogranular materials for 
weak tunneling conductance between the grains, gt < 1. We calculate the thermopower and figure 
of merit taking into account the shift of the chemical potential and the asymmetry of the density 
of states in the vicinity of the Fermi surface. We show that the weak coupling between the grains 
leads to a high thermopower and low thermal conductivity resulting in relatively high values of 
the figure of merit on the order of one. We estimate the temperature at which the figure of merit 
has its maximum value for two- and three-dimensional samples. Our results are applicable for 
many emerging materials, including artificially self-assembled nanoparticle arrays. 
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Introduction. — Thermoelectric materials with high 
efficiency are a major research area in applied physics 
dnd materials science for several decades now. Due 
to recent advances in nano-fabrication, these materials 
jiromise next generation devices for conversion of ther- 
mal energy to electrical energy and vice versa. Especially 
Suited for further improvement in efficiency are inhomoge- 
neous /granular thermoelectric materials [1] in which one 
dan directly control the system parameters. A measure 
for the performance or efficiency of a thermoelectric ma- 
terial is the dimensionless figure of merit, usually denoted 
as ZT, where T is the temperature. It depends on the 
thermopower or Seebeck coefficient, 5* and the electric, a 
and thermal, n conductivities, ZT S'^aT/n. [2,3] In 
Ref. [4] a figure of merit at 300i^ of 2.4 for a layered 
nanoscale structure and later in Ref. [5] a ZT of 3.2 at 
about 600i^ for a bulk material with nanoscale inclusions 
were reported. All these experimental achievements and 
technological prospects call for a comprehensive theory 
able to provide a quantitative description of thermoelec- 
tric properties of granular materials, which can in future 
serve as basis for a clever design of devices for a new gen- 
eration of nano-thermoelectrics. Recently we started this 
task by studying granular metals [6,7] and here we will an- 
swer the question to what extend weak coupling between 
grains is important in changing ZT . 

In this Letter we investigate the thermopower S and the 
figure of merit ZT of granular samples focusing on the case 
of weak coupling between the grains, gt < 1, see Fig. 1. 




Fig. 1: (Color online) Sketch of a nanogranular material show- 
ing typical electron (e) and hole (h) trajectories for asymmetric 
density of states in an electric field (E). The energy transport 
goes from the "hot" (H) to the "cold" (C) side. 



Each nanoscale cluster is characterized by (i) the charging 
energy Ec = e^ /{tro), where e is the electron charge, 
the sample dielectric constant, and a the granule size, and 
(ii) the mean energy level spacing 5. The charging energy 
associated with nanoscale grains can be as large as several 
thousands Kelvins thus the temperature range T < E'c in 
which our theory applies, is technologically relevant. 
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Model and Results. — The internal conductance of a 
grain is taken much larger than the inter-grain tunneling 
conductance, which is the standard condition for granu- 
larity. The tunneling conductance is the main parameter 
that controls macroscopic transport properties of the sam- 
ple [8]. In consideration of application to experiments [4,5] 
we restrict ourselves to the case where the tunneling con- 
ductance is smaller than the quantum conductance. 

The main results of our work are as follows, (i) We 
derive the expression for the thermopower S of granular 
materials in the limit of weak coupling between the grains 
taking into account a) the shift of the chemical poten- 
tial A/i ~ ai{v — \/2)To with v being the electron filling 
factor (\v\ < 1), fli a dimensionless numerical coefficient, 
and To = e'^ / {£r£,ioc) ~ is the characteristic tempera- 
ture scale with £^ioc — a/ h\{El /T'^gt) being the inelastic 
localization length [8] and b) the asymmetry of the density 
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Here d = 2, 3 is the dimensionality of a sample. We note 
that the thermopower S of granular materials is finite only 
if the shift of the chemical potential, A/i and the asymme- 
try of DOS, An are taken into account. Thus, for filling 
factor = 1/2 the thermopower vanishes. 

(ii) We obtain the expression for figure of merit of gran- 
ular materials in the limit of weak coupling between the 
grains as follows: 



ZT = 



2d^gte~V^ [Afi/T + AnTp]' 
7ig?(r/i?e)2 + j2Alph/ar-''{T/QDr-' ' 
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Here Iph is the phonon mean free path, O^) is the De- 
bey temperature, and the parameters A/i and An are de- 
fined above Eq. (1) [71 = 327r3/15, 72,20 = 12C(3)/7r, and 
72, 3D = Btt^/IS]. For two and three dimensional samples 
(d ~ 2, 3) the figure of merit ZT has its maximum value 
at temperatures T^ = To/[4:{d+ 1)^] [9]. We note that 

^2D ^ ^3D- 

Our main results [Eqs. (1) and (2)] are valid for tem- 
peratures T < Tq and weak coupling between the grains 
gt < 1. Under this conditions, the electric conductivity tr 
and thermal conductivity k = Kg + Kph, with Kg and Kph 
being the electron and phonon contributions to thermal 
conductivity, are given by expressions [8, 10-12] 
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We note that the thermopower S of lightly doped semi- 
conductors was investigated in the past [13-17]. How- 
ever, all previous studies were concentrated on the Mott 
variable range hopping (VRH) regime, with conductivity 
being a{T) - exp[-(TA//T)i/4], where Tm is the Mott 
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Fig. 2: (Color online) Result of simulations of the density of 
state n(e) vs. energy for different filling factors v and different 
system sizes. The simulations were done for a 2D Coulomb 
glass model at zero temperature [19, 21, 22]. We chose two 
large system sizes (A'^ = 500'^ and N = 800^) which show no 
size dependence for the shift of the chemical potential A/i and 
the asymmetry of the density of states An. Using simulation 
data we extract the numerical coefficients a\ and 02 in Eq. (1). 
We use the Coulomb glass model to simulate the whole weakly 
coupled granular sample. 

temperature [18]. In granular materials the Mott VRH 
regime is hard to observe. Indeed, the Efros-Shklovskii 
(ES) law [19,20] may turn into the Mott behavior with an 
increase of temperature. This happens when the typical 
electron energy e involved in a hopping process becomes 
larger than the width of the Coulomb gap Ac, i.e. when 
it falls into the flat region of the density of states where 
Mott behavior is expected. To estimate the width of the 
Coulomb gap Ac, one compares the ES expression for the 
density of states n(Ac) ~ (e^/e^)"'! Ad'*"^ with the DOS 
in the absences of the long-range part of the Coulomb 
interaction, ng. Using the condition n(Ac) ~ no we ob- 
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Inserting the value for the bare 



tain Ac 

DOS, no = l/EcSf", into the last expression we finally ob- 
tain Ac ~ Ec. This equation means that there is no flat 
region in the low temperature DOS. 

Thermopower. — Here we show that the ther- 
mopower S of granular materials is finite only if the shift 
of the chemical potential and the asymmetry of DOS are 
taken into account. In the following we outline the deriva- 
tion of Eqs. (1), (2). 

To calculate the thermopower of granular materials in 
the regime of weak coupling between the grains it is neces- 
sary to take into account electrons and holes because the 
contributions of electrons and holes cancel in the leading 
order. In general the thermopower is proportional to the 
average energy transferred by charge carriers and can be 
written as [13-17] 



5 = -^[(e-M>e + (e-M>J 



(4) 



Here the subscripts e and h refer to electrons and holes 
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Fig. 3: (Color online) Plots of figure of merit ZT vs tempera- 
ture T for two and three dimensional samples. Typical material 
parameters are given in the text. 

and /i = /i + A/i is the shifted chemical potential with 
A/1 = a\(y — l/2)ro being the chemical potential shift. 
The expression in the square brackets of the r.h.s. of 
Eq. (4) describes the average energy transferred by charge 
carriers (electron or hole) measured with respect to the 
shifted chemical potential /i. The average energy in Eq. (4) 
can be calculated as follows 

J de {e — jl)n{e — fl)f{e — jl)e 
J de n{e — )i)f{e — //)e 



Here /(e — p.) is the Fermi distribution function, A = 
\/TqT is the typical transfer energy in one hop, and n[e — 
jl) is the DOS. 

Equation (5) differs in several important ways from pre- 
vious works [13-17]: all considerations in the past were 
concentrated on the Mott VRH regime with filling fac- 
tor V = 1/2. In our consideration the average energy, 
{e — fi)e + {s — jj)h is zero for half filling factor. To obtain 
the finite result in Eq. (4) it is crucial to take into account 
the asymmetry of the DOS and the shift of the chemical 
potential. 

The DOS, n{e - fi) in Eq. (5) for /i - Ac < e < /i has 
the form (see [14, 15]) 

n(e) cx |e-/2|''-^[l- (£-/i)A?7,], (6) 

and is constant, tiq outside the Coulomb gap region, e < 
p, — Ac and s > fl, where Ac ~ Ec is the width of the 
Coulomb gap. The shift of the chemical potential A/i = 
/i — /i and the asymmetry of the DOS, An, are explicitly 
defined above Eq. (1). 

Simulations. — To support our choice for the expres- 
sion of the DOS 7i(e— //) in Eq. (6) we numerically compute 
the DOS for a 2D Coulomb glass model ^ (see Refs. [21,22] 

^For the 3D case we expect similar values for the numerical pref- 
actors ai^2i it was shown in [22] for the hard gap in the weakly 
disordered case. 



for details) at arbitrary filling factor using first princi- 
ples. The result of the simulations is shown in Fig. 2. 
These simulations clearly indicate that for a filling fac- 
tor v ^ 1/2, the DOS is asymmetric and the chemical 
potential is shifted. By fitting the obtained DOS by ex- 
pression (6) for different ly and linear regression of A/i(j/) 
and An{i>), we can identify the dimensionless coefficients 
oi ~ 2.9 and 02 ~ 0.4. We note that ai ^ 02 and there- 
fore the contribution to the thermopower caused by A/i in 
Eq. (1) is dominant. 

Now, we can calculate Eq. (5) and the analog contri- 
bution for holes using Eq. (6) for the DOS together with 
the numerically found values for ai_2- Finally one obtains 
the result for thermopower S in Eq. (1). Substituting 
Eqs. (1) and (3) into the expression for the figure of merit 
ZT = S'^aT/{Ke + Kph) one obtains Eq. (2). 

Discussions. — In Fig. 3 we plotted the figure of 
merit ZT for a two- and three-dimensional system, us- 
ing typical parameters for granular materials: grain size 
a ~ lOnm and phonon mean free path Iph ~ a ~ 
Inm. The latter is in general temperature dependent, 
which can be neglected for the temperatures and grain 
sizes under consideration. With Qjj ~ 450K, Tq ~ 
4000K, and gt ~ 0.1 at temperature T = 200K us- 
ing Eq. (2) we can simplify ZT ~ {2/j2,d){o-id)'^ gt{i^ — 
l/2fe-'/^{QD/TY-^{To/Tf. Using {u-l/2) = 0.2 
we estimate {ZT)2d ^ 0.60 and (ZT)3£, ~ 0.25 (see dots 
in Fig. 3). An important difference between 2D and 3D is 
that ZT is proportional to Iph in 3D, whereas it is inde- 
pendent in 2D. For the chosen value of Iph both ZT curves 
are of the same order. However, for larger Iph, ZT also be- 
comes larger in 3D which is related to the fact that the 
scattering of phonons in 3D is more effective than in 2D 
due to the larger number of grain boundaries. This leads 
to the suppression of Kph and to higher values of ZT. 

For hybrid structures of different types of grains this 
effect might be enhanced and will be subject of a forth- 
coming work. Last we mention that in the regime of VRH 
the contribution due to phonon drag to thermopower is 
absent [14,15]. 

Conclusions. — In conclusion, we studied thermo- 
electric properties of inhomogencous / nanogranular mate- 
rials at weak tunneling conductance between the grains, 
gt < 1- We calculated the thermopower S and figure of 
merit ZT taking into account the shift of the chemical po- 
tential and asymmetry of the density of states. We showed 
that the weak coupling between the grains leads to high 
thermopower and low thermal conductivity resulting in 
relatively high values of figure of merit. We calculated 
the temperature T* at which the figure of merit has its 
maximum value. 
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